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ABSTRACT 

Due to increasing^nced to account for the uncertainties in material properties, loading conditions, geometry'; 
-etc* T a methodology lias been developed to determine structural reliability aiui to assess the risk associated with 1 
it. The methodology consists of a probabilistic structural analysis by a probabilistic Finite element computer code 
'N-ESSUS (Nonlinear Evaluation of Stochastic Structures Under Stress)*and a generic probabilistic material 
properties model. The methodology is versatile and is equally applicable to high and cryogenic temperature 
structures. Results obtained demostrate that the whole issue of structural reliability and risk can be formally 
evaluated using the methodology developed which is inclusive of uncertainties in material properties, structural 
parameters and loading conditions. The methodology is described in some detail with illustrative examples. 



INTRODUCTION 


The probabilistic structural analysis method (PSAM) 1 has been developed to analyze the effects of fluctu- 
ating loads, variable material properties, and uncertain ties in analytical models, especially for high performance 
structures such as SSME turbopump blades. In the deterministic approach, uncertainties in the responses arc., 
not quantified and the actual safety margin remains unknown. Risk is calculated after extensive service experi- 
ence* However, probabilistic structural analysis provides a rational alternative method to quantify uncertainties 
in structural performance and durability. NESSUS is a probabilistic structural analysis computer code which 
integrates finite element methods and reliability algorithms, 2 * 3 capable of predicting the scatters of structural 
response variables such as stress, displacement, natural frequencies, buckling loads, etc. These are subsequently 
compared with their probable failure modes to assess the risk pf component fracture. Probable failure modes 
are included for different structures and their respective service environments. For example, failure events such 
as stress greater than strength, displacements exceed maximum allowables or avoidance of resonance are often 
used for the reliability assessment. Probability of occurrence of those failure events can be determined once the 
probability distributions of the requisite structural response variables are calculated by NESSUS. 

In undertaking a reliability/risk analysis, all suspected sources in uncertainties must be taken into account in 
order to control the probability of failure in service environments within an acceptable range. Structural reli- 
ability and risk obtained by a formal probabilistic methodology can be useful in evaluating the traditional design, 
setting quality control requirements, inspection intervals and retirement for cause. It can also be used to identify 
candidate material and design concepts in the absence of a technology base. In this report, the methodology 
developed to assess the structural reliability and risk/cost is described. 

CONCEPT OF PROBABILISTIC STRUCTURAL ANALYSIS 

In a probabilistic structural analysis, the primitive variables which define the structure have to be identified. 
These include temperature, material properties, structural geometry, loading conditions etc. and should be de- 
scribed by their respective probability distributions. A structural analysis performed by NESSUS with the pre- 
determined probability distributions of ail the primitive (random) variables will produce corresponding scatter 
(uncertainties) in the structural responses such as displacement, stress, natural frequencies etc. The concept is 
illustrated by Figure (1) where the structural model synthesizes the input uncertainties. 
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STRUCTURE OF NESSUS 


NESSUS consists of 3 major modules; NESSUS/PRE, NESSUS/FEM, and NESSUS/FPI. 

NESSUS/PRE is a pre-processor used for the preparation of the statistical data needed to perform the 
probabilistic finite element analysis. It allows users to describe the uncertainties in the structural parameters 
(random variables) at the nodal points of a finite element mesh. The uncertainties in these parameters are spec- 
ified over this mesh by defining the mean value and standard deviation of the random variable at each point, 
together with an appropriate form of correlation. Correlated random variables arc decomposed into a set of un- 
correlatcd vectors by a modal analysis. For strong correlation problems, the number of dominant random vari- 
ables in the set of uncorrelated vectors will be much less than that of correlated random variables. The 
computational time required for the analysis will also be reduced significantly. 

NESSUS/FEM is a finite element code used for the structural analysis and parameter sensitivity evaluation. 
It generates a database containing all the response information corresponding to a small variation of each inde- 
pendent random variable. The algorithm used in NESSUS/FPI requires an explicit response function in order 
to perform a reliability analysis. In complicated structural analysis problems, response can only be available 
implicitly through a finite element model. To overcome .this difficulty, the response function is expressed 
parametrically with this database. 

NESSUS/FPI (Fast Probability Integrator) is an advance reliability method 3 This module extracts the da- 
tabase generated by NESSUS/FEM to developed a response or a performance model in terms of uncorrelated 
random variables. The probabilistic structural response is calculated from the performance model. For a given 
response value, the probability of exceedence at this value is estimated by a reliability method, which treats the 
problem as a constrained minimization. This step is called a point probability estimation. The cumulative dis- 
tribution function is generated by running FPI at several response values. One alternative for generating the 
distribution function for any given response is to conduct a direct Monte Carlo simulation study. However, in 
general, it is very costly. NESSUS/FPI provides a method which not only produces a reliable distribution, but 
also requires less computing time than that for Monte Carlo simulation, especially in low probability regions. 

PROBABILISTIC MATERIAL PROPERTIES MODEL 

A generic material behavior model 4 is used to evaluate the scatters of material properties for the structures 
subjected to high temperature environments and high cycle loading conditions. The fundamental assumption for 
this model is that the material properties behavior can be simulated by primitive (random) variables. The general 
form of this model is shown in TABLE (1). The exponents n, p, q are determined from available experiment data 
or can be estimated from the anticipated material behavior due to the particular primitive random variable. 

APPLICATION TO A SPACE SHUTTLE MAIN ENGINE (SSME) TURBINE BLADE 

The methodology described previously will be applied to SSME turbine blade as shown in Figure (2). This 
blade is subjected to complex mechanical and thermal loads 5 ; centrifugal force, thermal loads and the differ- 
ential pressure across the airfoil. Centrifugal force is induced by rotational speed. Since it is difficult to maintain 
a constant rotational speed, the centrifugal force has to be considered as a random variable. Random thermal 
loads are due to combustion irregularities which causes a random temperature distribution in the blade. Differ- 
ential pressure is also random because of pressure fluctuation. Uncertainties in the blade geometry arise during 
the manufacturing process. The scatter in material properties is caused by nonuniformities in the material. 

The turbine blade is modelled by 40 4-node shell elements with 55 nodal points. In this study, 7 random 
fields are considered as listed in TABLE (2). In previous analysis, 6 the probability density function of material 
properties such as Young's modulus, thermal expansion coefficient and material strength are assumed. In the 
present study, they are simulated by the probabilistic material property model as defined in the TABLE (1). The 
material property is a function of its reference value, temperature, and fatigue cycles. The statistics of the prim- 
itive random variables are listed in TABLE (3). Since the material is a function of stress and stress is a function 
of material property in an implicit way, an iterative procedure is necessary to obtain a convergent solution for 
stress and material properties. The procedure is illustrated in TABLE (4). Young's modulus, thermal expansion 
coefficient and material strength are considered to be random variables. During the iteration process, the joint 
cumulative distribution functions of nodal stresses are also calculated in order to apply the generic probabilistic 
material property model on the entire blade. It is found that only a few iterations are needed for a convergence. 
The blade is also assumed to be subjected to 100,000 constant amplitude load cycles which degrade the modulus 
and strength. The critical points where the large displacement or high stress located are depicted in Figure (3). 
As shown in Figure (4), at the, root of the leading edge, the probability distribution of modulus is reduced sig- 
nificantly after the application of cyclic loads. Changes in the probability density function of tip displacement 
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with and without the application of cyclic loads are also shown in Figure (5). As expected, the tip displacements 
arc increased after the application of cyclic loads because the blade becomes softer. Material strength degrada- 
tion after cyclic loads is calculated with the converged stress by the probabilistic material property model as 
shown in Figure (6). Varying the number of cycles and repeating the procedure described previously, stress 
strength relationships are determined for a given cycle. These relationships are used to define a failure mode as 
the event when the stress is greater than the corresponding material strength. The probability of failure is then 
calculated by 


p f - ds)f a {x) (lx (i) 

where f a is the probability density function (pdf) of effective stress calculated by the probabilistic structural 
analysis using NESSUS. f s is the pdf of material strength simulated by the probabilistic material property 
model. From this analysis, a Risk-Fatigue cycle curve is developed for critical locations as shown in Figure (7). 
This curve is useful for assessing the risk of structural fracture. For instance, at a given acceptable risk level, the 
number of fatigue cycles to initial local failure can be determined. With this information available, criteria can 
be set for quality control, inspection intervals and retirement for cause. 

The risk assessment includes the initial cost and the consequential cost. The initial cost is defined as the cost 
for component service readiness which can be a function of several key design variables. The consequential cost 
is the cost incurred due to failure. Total cost is the sum of initial cost and a fraction of consequential cost as 
defined in Equation (2). The fraction is weighted by the probability of failure. 


C t = C; + Cf ( 2 ) 

where C f represented initial cost and Cf is the consequential cost if failure occurs. Since the lower initial cost 
is often associated with higher risk for the structural failure and higher initial cost will normally reduce the risk, 
the total cost can be minimized for an acceptable structural reliability. For example in Figure (8), the initial cost 
is a function of the quality of the material; more cost is needed for improved quality, yet the cost due to failure 
is reduced. There exists a point of diminishing return. This can be identified as the lowest total cost but with no 
loss of the structural reliability. A similar result can be obtained for the case of mean strength improvement as 
shown in Figure (9). 


SUMMARY 

* 

In summary, a reliability/risk cost methodology has been developed. It consists of a probabilistic structural 
analysis by NESSUS and a generic probabilistic material model. The methodology is versatile and equally ap- 
plicable to hot and cold structures where data is difficult to obtain. The methodology is demonstrated by using 
it to assess the risk associated with fatigue cycles to initiate local failure in SSME blades. 
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FIGURE (2) UNCERTAINTIES IN THE PROBABILISTIC STRUCTURAL ANALYSIS 
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FIGURE (3) SSME BLADE SHOWING LOCATIONS WHERE PROBABILISTIC 
STRUCTURAL RESPONSE WAS EYALUATIED 



(MEAN STRESS STRENGTH RATIO - 0.83) (MEAN STRESS STRENGTH RATIO - 0.68) 


PROBABILITY 
OP OCCURRENCE 



FIGURE (4) PROBABILISTIC MODULUS SIMULATED BY USING THE GENERIC 
PROBABILISTIC MATERIAL PROPERTY MODEL ( MPSI ) 


PROBABILITY 
OF OCCURRENCE 



FIGURE (0) PROBABILISTIC DISPLACEMENTS CALCULATED BY NBSSUS (10 3 INCH) 


(MEAN STRESS STRENGTH RATIO « 0.63) 


(MEAN STRESS STRENGTH RATIO - 0.68) 


PROBABILITY 
OF OCCURRENCE 



FIGURE (6) PROBABILISTIC FATIGUE STRESS/STRENGTH SIMULATED BY USING 
THE GENERIC PROBABILISTIC MATERIAL PROPERTY MODEL ( KSI ) 


v 


159 



PROBABILITY 
OP FAILURE 


FIGURE (7) 
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FIGURE (9) TOTAL COST AND STRUCTURAL RELIABILITY QUANTIFIED 
IN TERMS OF QUALITY CONTROL (GIVEN MEAN STRENGTH) 
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PRIMITIVE VARIABLES 
M p = MATERIAL PROPERTY 
T = TEMPERATURE 
S = STRENGTH 
a = STRESS 

N m = MECHANICAL CYCLES 


SUBSCRIPTS 

F = FINAL CHARACTERISTIC VALUE 
0 = REFERENCE PROPERTY 


TABLE(l) GENERIC PROBABILISTIC MATERIAL PROPERTY MODEL 
IN TERMS OF PRIMITIVE VARIABLES 


RANDOM P I LEDS 

NUHRKR OF 
DEPENDENT 
R.V. 

MEAN 

STANDARD 
DEVIATION 
(OR COV 3 

CORRELATION 

l.ENCTII 

NNMDER OF 
INDEPENDENT 
R.V . 

X COORDINATE 

55 

deterministic 

COORDINATE 

0.01 

l n 

5.0 

13 

Y COORDINATE 

55 

DETERMINISTIC 

COORDINATE 

0.01 

in 

5.0 

13 

Z COORDINATE 

55 

DETERMINISTIC 

COORDINATE 

0.01 

In 

5.0 

13 

TEMPERATURE 

55 

STEADY STATE 
TEMPERATURE 

60 

° F 

3.0 

22 

MODULUS 

AO 

23 HKS I 

0.10 

(COV) 

3.0 

16 

PRESSURE 

36 

STEADY state 
PRESSURE 

0.20 

(COV) 

n.o 

36 

ROTATION SPEED 

1 

A OOOO RPM 

0. 01 

(COV) 

N/A 

1 


TABLE (2) RANDOM INPUT DATA 


VARIABLE 

DISTRIBUTION 

MEAN 

STANDARD DEVIATION 


TYPE 


(VALUE) 

(% OF MEAN) 

r F 

NORMAL 

2750 *F 

61.4 ‘F 

2.0 

1 0 

NORMAL 

68 *F 

2.04 *F 

3.0 

s F 

NORMAL 

212.0 ksl 

10.6 ksl 

5.0 

°o 

CONSTANT 

0 

0 

0 

n mf 

LOGNORMAL 

10 8 

5X10® 

5.0 

n mo 

LOGNORMAL 

10 3 

60 

5.0 

n 

NORMAL 

0.25 


3.0 

P 

NORMAL 

0.25 


3.0 

Q 

NORMAL 

0.25 


3.0 


TABLE(3) PRIMITIVE VARIABLE PROBABILITY DISTRIBUTIONS FOR 
PROBABILISTIC MATERIAL PROPERTY MODEL 
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ITERATION 0 

STEP 1: CALCULATE THE PROBABILITY DISTRIBUTIONS OF 

MATERIAL PROPERTIES USING THE GENERIC 
MATERIAL MODEL WITHOUT STRESS AND FATIGUE 
CYCLES 

STEP 2: CALCULATE THE NODAL STRESS BY NESSUS 

ITERATION 1 

STEP 1: CALCULATE THE PROBABILITY DISTRIBUTIONS OF 

MATERIAL PROPERTIES USING COMPLETE GENERIC 
MATERIAL MODEL 

STEP 2: CALCULATE THE NODAL STRESS BY NESSUS WITH 

UPDATED MATERIAL PROPERTIES 

ITERATION 2 

(REPEAT THE ITERATIONS UNTIL THE PROBABILITY 
DISTRIBUTIONS OF MATERIAL PROPERTIES AND 
STRESS HAVE CONVERGED) 


TABLE (4) PROBABILISTIC STRUCTURAL ANALYSIS BY NESSUS 
& A GENERIC PROBABILISTIC MATERIAL MODEL 



